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Abstract 

Numerical simulations of a model of plane Couette flow focusing on 
its in-plane spatio-temporal properties are used to study the dynamics 
of turbulent spots. While the core of a spot is filled with small scale 
velocity fluctuations, a large scale flow extending far away and occu- 
pying the full gap between the driving plates is revealed upon filtering 
out small scales. It is characterized by streamwise inflow towards 
the spot and spanwise outflow from the spot, giving it a quadrupolar 
shape. A correction to the base flow is present within the spot in the 
form of a spanwise vortex with vorticity opposite in sign to that of 
the base flow. The Reynolds stresses are shown to be at the origin of 
this recirculation, whereas the quadrupolar shape of the in-plane flow 
results from the transport of this recirculation by the base flow that 
pumps it towards the spot in the streamwise direction and flushes it in 
the spanwise direction to insure mass conservation. These results shed 
light on earlier observations in plane Couette flow or other wall flows 
experiencing a direct transition to turbulence by spot nucleation. 

1 Introduction 

Being stable against infinitesimal perturbations for all Reynolds numbers, 
plane Couette flow (pCf), the shear flow between two parallel plates mov- 
ing in opposite directions with velocity ±f/p, is the prototype of flows that 
require localized flnite amplitude disturbances to be pushed towards a tur- 
bulent regime. The transition is thus characterized by the nucleation and 
nonlinear growth of domains of turbulent flow, separated from laminar flow 
by sharp fronts and called turbulent spots (e.g., [1-5]). This kind of transition 
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is not restricted to pCf but is also present in plane Blasius (boundary layer) 
flow [6,7] or plane Poiseuille flow [8]. A review of some relevant laboratory ex- 
periments is given by Henningson et al. [9] and of their numerical counterpart 
given by Mathew & Das [10]. In practice, direct transition to turbulence via 
spots can be expected whenever no low-Reynolds number instability of iner- 
tial origin exists, whereas turbulent solutions to the Navier-Stokes equations 
may exist and compete with the laminar base flow at moderate Reynolds 
number [11, Chap. 6, §6.3]. 

Growing turbulent spots in pCf have been studied both experimentally [1- 
5] and numerically [12-14]. In their pioneering direct simulations of Navier- 
Stokes equations with realistic no-slip boundary conditions, Lundbladh & 
Johansson [12] pointed out that (i) the wall-normal velocity component — 
typical of internal irregular small scale structures — faded away outside the 
spot but (ii) slowly varying in-plane velocity components extended far outside 
with an inwards streamwise motion towards the spot at the streamwise edges 
and an outward spanwise motion at its spanwise edges. These observations 
were made by low-pass Gaussian flltering the small scales of the velocity 
fleld at mid-gap. Tillmark [5] conflrmed them experimentally by detecting 
the outwards spanwise component that developed over the full gap between 
the plates. 

More recently, Schumacher & Eckhardt [14] re-investigated the growth 
of turbulent spots by means of direct numerical simulations but using un- 
realistic free-slip boundary conditions at the plates. By averaging the flow 
field between the two plates, they also observed that the turbulent spot was 
accompanied by an overall spanwise outflow and streamwise inflow, which 
they termed quadrupolar . 

Spots seem to behave as obstacles in the base flow [3,7, 15]. Accord- 
ingly, they introduce additional pressure flelds induced by the distribution 
of Reynolds stresses associated with the small scale fluctuations inside the 
spot and generating the large scale flows. A similar interpretation was put 
forward by Hayot & Pomeau [16] who introduced a back-flow to explain the 
organization of spiral turbulence in cylindrical Couette flow [17], with possi- 
ble application to the banded turbulent regime discovered more recently in 
pCf [18] and numerically studied by Barkley & Tuckerman [19]. 

Previous experimental studies by Bottin et al. [20] have shown that, in 
the lowest part of the transitional Reynolds number range, flow patterns of 
interest extend over the full gap. We take advantage of this observation 
to study the dynamics of spots using numerical simulations of a previously 
derived model of pCf shown to display sufficiently good properties for this 
purpose [21]. The model is sketched in ^ and completed in the Appendix. 
Typical results of simulations are presented in ^ emphasizing the output of 
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the filtering procedure: (i) the in-plane quadrupolar fiow outside the spot 
and (ii) a spanwise recirculation cell inside. These observations are then in- 
terpreted in ^ where the generation of these two large scale fiow components 
is explained in terms of Reynolds stresses averaged over the surface of the 
spot. In the concluding section, we summarize our results and point to their 
relevance to the interpretation of previous observations in other wall fiows of 
less academical interest, such as plane Poiseuille [23] or Blasius fiows [24]. 



The model used here is an extension to realistic no-slip boundary condi- 
tions of an earlier model proposed by one of us [22] for unrealistic free-slip 
boundary conditions. It is derived in [21] from the Navier-Stokes equations 
through a systematic Galerkin method involving expansions in terms of poly- 
nomials, functions of the cross-stream coordinate y multiplied by amplitudes 
describing the in-plane (x, z) space dependence of the full velocity field. The 
equations are written for the perturbation {u' , v' , w' , p') to the laminar basic 
fiow t/bX, where x denotes the streamwise direction, i.e. u = Uh{y) + u'; v' 
and w' denote the perturbations in the cross-stream and spanwise directions, 
respectively, p' being the pressure perturbation. Lengths are scaled by the 
half-gap between the plates h, and velocities by Up so that the time scale is 
h/Up. The control parameter is the Reynolds number defined as R = Uph/u, 
where u is the fiuid's kinematic viscosity, and the dimensionless base fiow 
profile reads f/b(?/) = y foT y E [—1, 1]. 

In accordance with experimental observations [20], truncation of the 
Galerkin expansion at lowest consistent order is performed, reducing the 
set of basis functions to: 



u'{x,z,t,y) = Uoix,z,t)B{l-y') + Ui{x,z,t)Cyil-y^), (1) 



w'{x,z,t,y) = Wo{x,z,t)Bil-y') + W,ix,z,t)Cy{l-y'), (3) 



where A, B, and C are normalisation constants. These expressions are in- 
serted in the continuity and Navier-Stokes equations, and projections of 
the results on the same basis functions using the canonical scalar product 
(/; 9) = I-i f{y)9{y) d?/, are performed, which yields a set of coupled partial 
differential equations. For example, the projection of the continuity equation 



only display the equation for the amplitude Uq of the streamwise velocity 
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component which is even in y. 

dtUo + Nu, = -d^Po - aAUi - a^Vx + i?"' (A - 70) , (5) 

where A = dxx + dzz and with: 

7V^„ = ai(t/o5,f/o + IVo5.f/o) + ^«2(f/i5.f/i + + W^d^W^) , (6) 

just to show that each equation has the form expected for a hydrodynamic 
problem. In particular, nonlinearities have the same structure as the classical 
advection term v ■ Vv. In the same way, the last term in ([5]), with the factor 
R~^, accounts for the viscous dissipation associated with the cross-stream 
parabolic (70) and in-plane dependencies of Uq. This flow component can 
straightforwardly be identified as the streamwise streak amplitude, so that 
the source term —a2Vi on the r.h.s. of ([5]) accounts for the lift-up mechanism 
since Vi is the cross-stream velocity fluctuation. The physical role of the 
linear term —aid^Ui will be considered later. 

On general grounds, the Reynolds-Orr equation governs the perturbation 
energy E{t) = | dV, where V is the volume of the domain. It 

can be symbolically written as = P—D, where P is the energy production 
issued from the interaction of the perturbation with the base flow Uh{y) = y, 
P = — Jyu'v'-^UhdV, and D is the dissipation due to viscous effects. In 
our model, one readily gets -P = — Xs xUqVi dS, where S is the surface of 
the domain and x is a positive constant. Since Vi generates Uq through the 
lift-up mechanism, regions where the Reynolds stress —UoVi is positive, thus 
destabilizing the base flow and contributing to the turbulence production, 
are those with Uo > and Vi < or the reverse, which obviously correspond 
to Q2 and Q4 events identified in the literature, see for example [25]. 

The main limitation of the model comes from its low order truncation. In 
fact, expressions ((iHHl) are only the first terms of series expansions and the 
derivation of models truncated at higher and higher orders remain possible. 
Up to now, this has not been done for several converging reasons, the main 
ones being that (z) the Reynolds number range we are interested in corre- 
sponds to the lower part of the pCf's transitional regime where departures 
from laminar flow are known to occupy the full gap [20], [ii) Ui already con- 
tains the lowest order non-trivial correction to the base flow thought to be 
important in the discussion of the laminar-turbulent coexistence [16]. Ac- 
cordingly we believe that the lowest order model is sufficient to account for 
the large scale features present in the experiment at least at a qualitative 
level, the alternative being to turn to direct numerical simulations and not 
to consider a much more cumbersome higher order model. The discussion in 
§4] supports the validity our approach a posteriori. 
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3 Numerical simulations of turbulent spots 



Our model was integrated on a rectangular {x, z) domain with periodic 
boundary conditions, while being written for stream-functions \E'o; and 
velocity potential $i related to the velocity amplitudes through: 

U, = Uo-d,^,, Wo = W, + dj^o, (7) 
f/i = f/i + 9,<l>i - , Wi = Wi + d,^i + d,^i, Vi = ^^i/(3. (8) 

A standard, Fourier based, pseudo-spectral code was implemented with non- 
linear terms and linear non-diagonal terms (e.g. —aidxUi — 02^1 in (j5])) 
evaluated in physical space and integrated in time using a second order 
Adams-Bashforth scheme. The necessary introduction of Uq,. . . is com- 
mented upon in the Appendix. Simulations were performed in a domain 
of size {Lx X Lz) = (128 x 128) with effective space steps 6x = 6z = 0.25 
and 6t = 0.01. These values were retained as a good compromise between 
accuracy and the possibihty to let sufficiently wide systems evolve over suf- 
ficiently large periods of time [21]. Concerning the accuracy problem, it 
should be noted that small-scale in-plane structures are pieces of streaks and 
streamwise vortices with typical size larger than 3, which makes more than 10 
collocation points per structure. Smaller time steps did not produce results 
different from those shown here during comparable time lengths. 

As an initial condition, we took localized expressions for "^q, \E'i, and $1: 

^o{x, z,t = 0) = ^i(x, z,t = 0) = $1(0;, z,t = 0) = A exp~(^'+^')/^ 

where A is an amplitude and S is the size of the germ. Parameters A = 5 
and S = 2 were found efficient in generating turbulent spots for R = 250, 
well beyond Rg ~ 173, above which sustained turbulence is expected in our 
model [21]. In practice, due to the highly unstable characteristics of the flow 
at such values of R, the apparent simplicity of the initial condition played 
no role after a few time units. 

Spots are best illustrated by their most spectacular feature, namely their 
streamwise streaky structure [1,2,7,8]. In turn, the latter is best visualized 
from the amplitude Uq since streamwise streaks are easily identified as regions 
where \Wo\ -C |^7o| alternating in the spanwise direction, and since Uq is 
associated to velocity perturbations that are maximum in the mid-gap plane 
y = 0. Figure [1] displays gray-level snapshots of Uq at different times after 
launching. Denoting by (xc, zq) the in-plane coordinates of the center of the 
spot we see that, contrasting with the cases of plane Poiseuille or boundary 
layer flows, the spot does not drift due to the absence of mean advection. 
One can also notice its overall ovoid shape with dominant negative values 
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Figure 1: Growth of a turbulent spot at i? = 250 in a wide domain (L^. xL^ = 
128 X 128). Field of amplitude Uo{x,z,t) in gray levels at t = 50, 150, 250 
and 350 (from left to right and top to bottom). The whole domain becomes 
uniformly turbulent at t ~ 700. 



(dark gray) for x > xc and positive values (light gray) for x < xq- Regions 
where Uq is positive correspond to high and low speed streaks for ?/ > and 
y < 0, respectively, which compares well with the experimental observations 
in [20]. 

In the sequel, we study the state at t = 150 but results and conclusions 
are identical at different times. The complete field (f/o, Wq) corresponding to 
this reference state is displayed in Figure O Except in the very center of the 
spot that looks rather messy, streamwise structures are easily recognized but 
the trace of the large scale quadrupolar flow, of main concern in the present 
paper, is already visible. 

As done by Lundbladh & Johansson [12], we now proceed to the elimina- 
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Figure 2: Streak flow field (f/o, 1^0) at t = 150. 



tion of small scales using a Gaussian filter in spectral space: 

t(A;,, K) = Z{h, h) exp[~{kl + kl)/{2af] , (9) 

where the hat denotes the Fourier transform of any quantity Z = Uo,- ■ ■ ■ 
In physical space, this corresponds to a convolution with a kernel oc 
exp {^—(y\JCx + Q) where a is the parameter controlling the width of the do- 
main over which the small scales are smoothed out by the operation. Small 
scales, indicated by superscript 's', are recovered afterwards from the relation 

= z -~z. 

The diameter of the Gaussian averaging window has to be chosen in 
accordance with the size of the modulations to be eliminated, here the small 
scale streaks with spanwise wavelengths of the order 3-6 as can be guessed 
from Figure O We used a = tt/II, but the results were found to be rather 
insensitive to this choice provided that a is sufficiently small. 

As seen in Figure 0, this filtering procedure yields a clear picture of the 
fiow outside the spot: the overall pattern formed by the in-plane components 
Uq and Wo has a quadrupolar aspect that could already be guessed from the 
consideration of the unfiltered stream-function whose Laplacian is related 
to its vortical contents. In what follows, we term drift flow the large-scale 
velocity field {Uq, Wq) with Poiseuille-like cross-stream profile by analogy 
with the case of Rayleigh-Benard convection where a fiow with the same 
global features was introduced by Siggia & Zippelius [26]. 

Figure m displays the velocity components associated to the fields \I'i, $1. 
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Figure 3: Top: level lines of averaged velocity components Uq (left) and Wq 
(right), illustrating large scale streamwise inflow and spanwise outflow around 
the spot. Bottom, left: representation of this flow as vectors. Bottom, right: 
level lines of the unfiltered stream function \E'o. 



The distribution of the amplitude of Vi, displayed in the left panel, represents 
an average wall- normal motion which is maximum in the mid-plane y = 0, 
positive on the right of the spot's center x > xc and negative on its left. In 
turn, the flow (f/i, Wi) shown in the right panel consists in a region centered 
around the spot where \Ui\ ^ \Wi\ and Ui < 0. This structure is easily 
interpreted as a wide spanwise recirculation cell with vorticity opposite in sign 
to that of the base flow. It is further reminiscent of what can be deduced 
from DNS simulations of Lundbladh and Johansson [12], as displayed in seen 
their Fig. 9. 

In Figure [5] (a) we display the profiles of Uo and U i along a streamwise 
line going through the center of the spot. The dashed line corresponds to Uq 
and clearly points out the inwards character of the drift flow. In contrast. 
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Ui (solid line) presents a deep trough at the location of the spot. At the 
spot's center where f/o — 0, the superposition of the perturbation u' = 
UiCyil — y"^) and the base flow f/b(l/) = Hi shown in Figure [5] (b), displays 
the characteristic S shape of the turbulent velocity profile expected for pCf. 
The presence of the spot thus locally increases the wall friction. At different 
positions inside the spot, where f/o 7^ (and Wq 7^ 0), the full superposition 
U{y) = y{l + UiC{l —y'^)) + UoB{l — y^) leads to asymmetric mean velocity 
profiles (Fig. [5](c) for point xl and (d) for point x-ei) that are reminiscent of 
the averaged profiles obtained by Barkley & Tuckerman in their simulations 
of the banded regime of turbulent pCf [19]. 



4 Generation of large scales from small scales 

The mechanism driving the quadrupolar drift flow is discussed in terms of 
equations obtained by filtering from the model's equations, as described in the 
Appendix. We focus on the slowly varying quantities Aq = A^^o, ^1 = A^'i, 
and A2 = A$i, driven by Bi = -(U^f where ^ = a2(/5 + P") > and 
B2 = ai{U^y - (W^y + a2iUfy - (Wiy. The latter quantities represent 
the components of the Reynolds stress tensor [27] which do not average to 
zero over the surface of the spot {Bi corresponds to the energy extracted from 
the laminar flow and B2 mostly to the energy contained in the streamwise 
streaks) . 

Introducing slow variables X and Z whose rate of change is inversely 
proportional to the width of the window that is dragged over the data upon 
averaging through ©, one can observe that, in the equations, the quantity 
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Figure 5: (a) t/i (red-solid) and (blue-dashed) as functions of coordinate 
X along the streamwise center-line, (b-d) Full average streamwise velocity 
profiles U{y) at x = xc (b), x = xl (c) and x = xr (d); the laminar profile 
f/b(?/) = y is indicated by a dashed-dotted line. 



B\ appears with one derivative in X or Z less than i?2, due to the fact that 
B\ substitutes one in-plane differentiation by a cross-stream OiV) differenti- 
ation. Further assuming that the spot is in a quasi-steady state [pt ~ 0) and 
that space derivatives are negligible when compared to OiV) constants when 
operating on the same quantities, at lowest significant order one can simplify 
Equations ( fT3] - [T^ to read: 

i?-So^ = «i (i^zAa - dxA^ , (10) 
R-'^iA, = dzB^-a^dxAo, (11) 
R-'liA2 = -dxBi. (12) 

The structure of this system invites one to examine the shape of the dominant 
Reynolds stress contribution Bi as a function of the slow variables. Figure [6] 
displays the averaged Reynolds stress field associated with the small scales 
—UqVi. As could be anticipated the latter is positive under the spot and one 
can furthermore observe its single-humped shape that, following Li & Widnall 
[15] who developed a similar approach for spots in plane Poiseuille flow, can 
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Figure 6: Distribution of the averaged Reynolds stress field —UqVi (left) and 
its variations along streamwise (solid line) and spanwise (dashed hne) cuts 
through the maximum of the distribution taken as the center of the spot at 
xc = 64, zq = 60 (right). 

be modelled as a Gaussian function of the form exp[— (X^ + Z^)/2]. This 
assumption will help us to make an educated guess about the mechanisms at 
work. 

Considering first Equation ffT^ . from the third equation in ([H]), i.e. 
Vi = A^i/P, we obtain that the contribution to Vi generated by Bi is 
~ Xexp[— (X^ + Z^)/2], i.e. a pattern with a positive hump for X > and 
a negative one for X < 0, resembling that in Figure H] (left). This velocity 
component forms with Ui a large scale recirculation loop. As seen from the 
first equation in ([H]), Ui contains two contributions of potential and rota- 
tional origins, respectively. In the neighborhood of the X axis, the variation 
of dxBi is dominated by its X dependence so that A2 = {dxx + dzz)^i ~ 

dxx^i = -dxBi and, accordingly, dx'^i Bi exp[-(X2 + Z'^)/2\. 

As to the rotational contribution —dz^i, from (ITTl) and forgetting the cou- 
pling with Aq (which is of higher order owing to the way it is generated from 
Ai and A2), we have similarly Ai = {dxx + dzz)^i ~ dzz^i ~ dzBi, hence 
—dz^^i ~ —Bi so that it adds constructively to the potential part. The 
resulting Ui closes the recirculation loop as inferred from Figure H] (right). 

Inserting Ai ~ dzBi and A2 ~ —dxBi in (ITOl) we obtain a right hand 
side in the form — XZexp[— (X^ + Z'^)/2] for A\I'o which is the vorticity 
contained in the (?7o, Wq) velocity field. This field displays four lobes with 
alternating signs. An approximation to the large scale drift flow along the 
axes can easily be obtained. Indeed, Uq can be obtained from Uq = —dz'^o 
by integrating Aq = {dxx + dzz)^o over Z and neglecting dxx^o since 
\i/o varies much less with X than with Z along the X axis. We obtain 
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f/o ~ — Xexp[— (X^ + 2'^)/2] which accounts for the observed inward flow 
along the streamwise center-line of the spot. The same argument can be 
transposed for the spanwise direction (now \t'o varies most rapidly in the X 
direction, which makes dzz^o negligible and eases the integration over X), 
yielding Wq ~ Z exp[— (X^ + Z^)/2] which similarly accounts for the outward 
flow along the spanwise center-line. Notice however that this solution is too 
approximate to fulfil the continuity condition accurately since computing 
dxUQ + dzWo leaves a residual of the form (X^ - Z2)exp[-(X2 + Z2)/2], 
though the main contribution in exp[— (X^ + Z'^)/2\ is nicely compensated 
near the origin where the Gaussian is at its maximum. At any rate the chosen 
shape is only a simplifying assumption. 

Physically, the spot is thus characterized by a mean correction to the base 
flow (represented in the model by f/i < 0) itself generated by a wall normal 
velocity component (here V^i) and forming a large recirculation loop. In turn, 
the transport of that mean correction (here UiCy{l—y'^)) by the base flow ap- 
pears to be a source term for the large scale drift flow (here {Uq, Wq)) whose 
pattern is enslaved to its streamwise gradient, balancing viscous forces and 
inertia (according to R~^'^qUq + aidJJi ~ 0) and expressing flow continuity 
(a,f/o + 9,1^0 = 0). 

5 Conclusion 

In this paper, we have studied the large scale structure of the flow inside and 
around a turbulent spot in a transitional pCf model focusing on the in-plane 
dependence of a small number of velocity amplitudes [21]. The approach 
is supported by the qualitative consistency between previous experimental 
results in the transitional regime [20] and our own numerical simulations of 
the model. 

Inside the spot, we flnd a wide spanwise recirculation loop with vorticity 
opposite in sign to that of the base flow. In particular, a patch of streamwise 
correction counteracting the base flow is observed, giving a S shape typical of 
turbulent flows to the velocity proflle inside the spot. A reduced model PH - 
\n\ links this recirculation to Reynolds stresses —UqVi generated by the small 
scale fluctuations. Outside the spot, the existence of an inward-streamwise 
outward-spanwise quadrupolar drift flow has been pointed out, the origin of 
which is attributed to a linear coupling with this recirculation and linked 
to linear momentum conservation through ffTUl) . By simply assuming that 
the region where the Reynolds stresses contribute to the turbulent energy 
production (i.e. —iU^Vi > 0) is one-humped with localized support, the 
main features of the large scale flow extracted from numerical simulations by 
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filtering are recovered. In tliis approach, we only focused on the generation of 
large scales by small scales but considered neither (i) the interactions between 
small scales themselves nor (ii) the feedback of large scales on small scales. 
Closure assumptions are clearly needed in order to have a self-consistent 
theory, and especially to explain the sustainment of turbulence within a spot, 
problem (i), and its spreading as time proceeds, problem (ii). 

Owing to the general character of the argument leading to their existence, 
one might also expect to find these large scale corrections in and around spots 
developing in transitional shear flows other than pCf for which they have al- 
ready been accounted for [5,12,14]. Evidence of their presence can indeed be 
obtained from Figure 12 reporting numerical work of Henningson & Kim [23] 
on plane Poiseuille flow and from Figures 6 and 9 describing the result of 
ensemble averaging of turbulent spots in boundary layer flow with slightly 
adverse pressure gradient in the laboratory experiments of Schroder & Kom- 
penhans [24]. Despite its limited cross-stream resolution, our modeling of 
transitional plane Couette flow has thus been shown to provide valuable ex- 
planations to previous observations, which might call for new laboratory ex- 
periments since, besides the theoretical challenge of understanding laminar- 
turbulent coexistence in detail, the problem of the transition to turbulence 
in wall flows has a great technical importance. 



A Model's equations and derivation of (1101- 

m 



As explained in the main text, the model is obtained by projecting the 
Navier-Stokes equations on the chosen basis (llHSl) with velocity perturba- 
tions expanded on the same basis. The set completing (jlj) and fl3|U|) reads: 



dtWo + Nwo = -d.Po - aid^Wi + R-\A - 70)^ 







Nwo = (yiW.Wo + Wod^Wo) + a2{UAW, + W^d^W, + (3'V,W,) , 

dtUi + Nu,= - d,Pi - a^^Mo + R'^i^ - li)Ui , 

Nu, = a2{UodMi + Uid^Uo + Wod^i + W,dM - P"ViUo) , 

dtW, + Nw,= - d,P, - a.d^Wo + R^' (A - 71) W, , 

Nw, = a2iUod^W, + U,d.,Wo + W^d^Wi + W^d^Wo - ^"ViWo) , 

dtVi + Nv, = -l3Pi + R-\A - ^[)Vu 

where A denotes the two-dimensional Laplacian dxx + dzz- Coefficients all 
derive from integrals of the form Jn,m = lo ^'"(1 -2/^)'"d?/ = EfeLo (m) 2fc+n+i • 
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3V15/14, as 



15/6, as 



We have: ai = 1/v^, 03 = ^27/28, ai = 
5x/T5/22, 70 = 5/2, 71 = 21/2, 7^ = /J^, = and /3" = \(3. 

The equations governing fields \l/o, ^l/i, $1, from which the velocity com- 
ponents derive through fl7|H|) . are obtained in the usual way by differentiating 
and cross- subtracting or adding the previous equations. They read: 



(9i-i?-i(A-7o))A^o 

= {d.Nu, - d^Nwo) + ai(|a,A$i - a.A^i) 

dt{A - P^) - R'\A^ - 2/3'A + ^^p^)] A$i 



P\d,Nu,+d,Nw,)-PAN. 



Vi ■ 



(13) 
(14) 

(15) 



The introduction of averaged quantities Uq, Wq, Ui, and Wi in ([7]) and ([H]) 
is forced by our choice of periodic boundary conditions, otherwise the pos- 
sibility of a uniform velocity correction corresponding to linearly increasing 
potential/stream functions would be overlooked. They are governed by: 
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a2(/3 + /5")fW-7i^"'t/i, 
a2{P + (3")WoVi--iiR-^Wi 

a2(/3-/3')tW-7oi?"'f/o, 
«2(/3-/3')W^i-7oi?"'W^o, 



where the wide tildes mean averaging over the whole domain. Among this set 
of equations, the first one is the most relevant since it precisely corresponds 
to the expected mean fiow correction. Quantity 02 (/3 + 13") was denoted ^ in 
the text. 

It was observed in Figure [T]that the fiow within the turbulent spot resem- 
bles developed turbulent fiow, see also [9,15]. Accordingly, one obtains that 
the only contributions to the averaged equations come from the terms that 
keep a constant sign over the surface of the spot, namely the main Reynolds 
stress term —UqVi associated with energy extraction from the mean fiow and 
the other terms Uq, Wq, f/f , and Wi. Equations (|T^1 - [T^ then reduce to: 



(9t-i?-i(A-7o))A^o 
{dt - R-\A - 7i))A^i = i^^UoV; 



9,A^ij, (16) 
(17) 

^^(A - 0') - R-'{A' - 20' A + 71/3^) A$i = 0'd, ( -^f/oVT) , (18) 



+ ai (|9,A<l>i 
- aid^A^o , 
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with ^ = (/5 + P")- Following Li & Widnall, we then split the velocity 
components into small and large scales, i.e. Uq ^ Uq + Uq, etc., and only 
keep the contribution to the Reynolds stresses coming from the small scales. 
This leads to the same set of equations as above except that Uo, Ui, ... are 
replaced by their small scale parts Uq, Uf, .... 
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